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We investigate the magnetic reorientation in a two-dimensional anisotropic antiferromagnet due 
to a transverse magnetic field. Using a many-body Green's function approach, we show that the 
magnetization component perpendicular to the applied field (and along the easy-axis of the anti- 
ferromagnet) initially increases with increasing field strength. We show that this unexpected result 
arises from the suppression of quantum and thermal fluctuations in the antiferromagnet. Above the 
Neel temperature, this effect leads to a reappearance of a magnetic moment along the easy-axis. 

PACS numbers: 75.70.Ak, 75.30.Ds, 75.50.Ee, 75.25.+Z 



I. INTRODUCTION 

The magnetic properties of two-dimensional (2D) an- 
tiferromagnets (AFM) have been extensively studied in 
the pasti2iil^ Examples of 2D AFMs that have been 
recently investigated include the manganites which ex- 
hibit a colossal magnetoresistance,^ the vanadates,^ and 
the undoped parent compounds of the high-temperature 
superconductorSiSi The latter are a prime example of 
weakly anisotropic 2D AFMs due to their small in-plane 
anisotropy and an even smaller interlayer coupling be- 
tween neighboring Cu02-planes. A finite anisotropy 
is necessary to stabilize the long-range magnetic order 
in 2D magnets at finite temperatures^ Even a rather 
weak anisotropy induces an ordering temperature of the 
same magnitude as the isotropic exchanged The prop- 
erties of the above materials have been intensely stud- 
ied theoretically within the framework of the Heisenberg 
modeliSiSi^iiSiiii^ii^ In particular, the magnetization of 
the isotropic 2D AFM as function of an applied magnetic 
field have been studied within spin-wave theory4*^ 

In this communication we study the properties of a 
2D anisotropic AFM with spin S' = 1/2 on a square lat- 
tice in a transverse magnetic field perpendicular to the 
easy axis of the anisotropy. To this end, we develop a 
many-body Green's function methodic that is based on 
the equation-of-motion formalism. Our results are two- 
fold. First, the staggered AFM magnetization along the 
easy axis increases with increasing strength of the trans- 
verse magnetic field. This effect is quite unexpected since 
the field is directed perpendicular to the magnetization 
component. We find, however, that the presence of a 
weak transverse field leads to the suppression of quan- 
tum and thermal fluctuations in the AFM, the former 
being responsible for the decrease of the zero tempera- 
ture sublattice magnetization from its saturation value, 
S = 1/2. While a similar behavior has been predicted 



for one-dimensional (ID) AFM spin chains,A2, its relation 
to our results is at present unclear due to the qualita- 
tively different nature of spin excitations in ID and 2D 
systems. Second, the AFM ordered spins take a non- 
coUinear canted orientation for any non-zero transverse 
field, inducing a non-zero magnetization component par- 
allel to the applied field. Note that this behavior is quali- 
tatively different from that of a 2D AFM in a longitudinal 
magnetic field parallel to the easy axis>ii^ In the latter, 
a canted spin configuration can only be reached via a 
phase transition into the so-called spin-flop-phase. 



II. THEORY 

We consider the anisotropic (XXZ-) Heisenberg Hamil- 
tonian 



n = 
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where is the spin operator with spin quantum number 
5=1/2 located on sites i of a square lattice, and J > is 
the isotropic exchange coupling between nearest neighbor 
(nn) spin pairs i and j. The easy-axis is modeled by an 
exchange anisotropy D > Q along the z-axis. We take 
the transverse magnetic field B — {B^, 0, 0) to be aligned 
along the x-axis, perpendicular to the easy axis. Similar 
results to the ones discussed below are expected for a 
single-ion anisotropy, although its physical origin differs 
from the exchange anisotropy considered in this study.— 
In order to compute the temperature- and field- 
dependence of the AFM's sublattice magnetization 
m.i{T, Bx), we employ a many-body Green's function 
approach, which is based on the equation-of-motion 
formalism>i& The non-coUinear magnetic structure which 
occurs in the presence of a non-zero transverse magnetic 
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field requires that two non- vanishing magnetization com- 
ponents have to be considered. This can be done either 
in the original spin coordinate systemji^ which however 
is analytically and numerically demanding due to the oc- 
curring zero-eigenvalue problem,^^ Therefore, we apply 
an approach identical to the one recently used to investi- 
gate the field-induced spin reorientation of 2D anisotropic 
ferromagnetS)S2i2i and to study the noncolUinear mag- 
netization of 2D isotropic antiferromagnets4i^ In this 
approach, the spins of each sublattice i are rotated lo- 
cally by angles 9i in such a way that in the rotated frame 
(primed spin operators) only a single non- vanishing com- 
ponent of the sublattice magnetization remains, i.e., 
TOi(r) = (Sf) / and (Sf) = {Sf) = 0. The sym- 
metry of the present case simplifies the calculation con- 
siderably by assuming equal magnitudes of the sublattice 
magnetization, ■m{T,Bx) — |mi(r, i?a;)| = \m.2{T, Bx)\, 
and canting angles 6*1 (T) = 9{T) and 6*2 (T) = tt - 9{T) 
with respect to the easy axis. 

Specifically, we consider the following commutator 
Green's functions in energy space, 



where 



(2) 



which we compute by the conventional equation-of- 
motion approach. We approximate higher-order Green's 
functions by using the Tyablikov decoupling for i ^ fc,^^ 



(5n((5f;5r'))^ = m,(T) G±.(c.) 



^ (3) 

For 2D fcrromagnets with a small anisotropy it has been 
shown that the Tyablikov decoupling (or random-phase 
approximation, RPA) yields almost quantitative results 
specifically for the magnetization and susceptibilities,^^ 
whereas the resulting free energy and the specific heat 
are less well described]^ For systems with spin quantum 
number 5 = 1/2 the magnetization can now be obtained 
from 



(4) 



a = Bx sme{T) 

+ q m(T) [J + D-{2J + D) sin^ e{T)] , (6) 

2 , 



b{k) = -m(T) 7(k) (J + D/2) sin^ 0(T) , (7) 
c(k) = m(T) 7(k) [J- (J + D/2) sin2 0(r)] , (8) 
rL(k) = [a±6(k)]2-c2(k). (9) 



7(k) — 2 {cos kx + cos ky) . 



(10) 



Here, q = 4 is the number of nearest neighbor sites, and 
the lattice constant oq is set to unity. The excitations of 
the system are represented by two branches of spin waves 
whose dispersions are given by ei.2(k). For B = 0, = 
the Neel temperature T/v is given by 
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(11) 



with a = q {J + D). Note that the "anomalous" Green's 
function {uj) appears due to the canted nature of 
the spin configuration in a transverse field. Its ap- 
pearance implies that in addition to the spin-flip term 
S~ S'^ also a term of the form S~ S~ is present in the 
HamiltonianiS^ As a consequence, the spin precession 
around the equilibrium direction is no longer spherical 
but elliptical. Also, the consideration of G~j{uj) guaran- 
tees that the Mermin- Wagner theorem is fulfilled for an 
easy-plane magnet (D < 0). 

In order to obtain the individual components of the 
magnetization, we need to compute the magnetization 
angle 9q{T, Bx)^ The latter can be obtained from the 
observation that in the rotated frame in equilibrium no 
torque is exerted on the magnetization (Sf), i.e., S^' 
is a constant of the motion, and d{Sf)/dt = 0. This 
requirement, in turn, is equivalent to the condition that 
the Green's function Gfj{oj) = {{Sf;S~ ))^ vanishes, as 
has been shown for the field-induced spin reorientation 
of a ferromagnetic monolayer^ Computing this Green's 
function by using the same treatment as for the {uj)- 
Green's functions, i.e., by application of the Tyablikov 
decoupling, one finds 



where N is the number of lattice sites, and the mo- 
mentum sum runs over the full Brillouin zone. The 
equal-time correlation function (S'j" )(k) is obtained 
after Fourier transformation into momentum space from 
Eq.(|2Jl via the spectral theoremii& Note that here the 
indices i, j = 1,2 refer to the two sublattices. We obtain 



(5r'5+')(k) = 



(T) 
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/3£2(k) 



coth 
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(5) 



sm9oiT,Bx) = 
Bx 



(12) 



q m{T,Bx){2J + D) ' Bx<qmiT,Bx){2J + D), 
1 , otherwise . 

This method to determine the angle 9o{T,Bx) was suc- 
cessfully applied by Schwieger et al.^^ and by Pini et al.^ 
for the case of the field-induced spin reorientation of a 2D 
ferromagnet. 

For r = we have compared this approach with the 
one where 9o{0,Bx) is determined from the minimum of 
the internal energy E{0, 9) = {H), and found satisfactory 
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FIG. 1: (a) Magnetization m{T,Bx) and (b) equilibrium an- 
gle 6o{T,Bx) as functions of the transversal field Bx for dif- 
ferent temperatures T below and above the Neel temperature 
Tjv. The temperatures and interactions are given in units of 
the isotropic exchange J. We have assumed a spin quantum 
number S — 1/2 and an anisotropy D/J = 0.01. 



agreement. In fact, the relation Ea. H12|) for 6o{T,Bx) 
is identical to the one calculated from the free energy 
-FmfaCT, 0) within a single-site mean field approximation. 
This approximation is easily obtained within our theoret- 
ical approach by neglecting the spin- flip term S^SJ in 
Eq.lP, or putting 7(k) = in Eqs.Q-® 



III. RESULTS 

Unless stated otherwise, we consider in the following a 
weak anisotropy D/J — 0.01, yielding a sublattice mag- 
netization of m(r = 0,B^ = 0) = 0.381 at T = 0. A 
magnetic anisotropy is necessary since for an isotropic 
system (i.e., for D ^ = 0), our approach satisfies 
the Mermin-Wagner-theorem.^*'-^^ and we have m = 
at T 7^ 0. For the isotropic 2D AFM (D = 0) on 
a square lattice we obtain m{T = 0, = 0) ~ 0.36, 
which compares reasonably well with the commonly ac- 
cepted value of 0.307.^'* The latter value is also obtained 
within the Holstein-Primakoff-approximation. The Neel 
temperature Tn for D/J — 0.01 is calculated to be 
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FIG. 2: Staggered magnetization component \mz{T,Bx)\ 
along the easy axis as function of the transversal field Bx 
for different temperatures T. For the other denotations we 
refer to Fig. 



Tn/J — 0.472, while a quantum Monte Carlo calculation 
yields a larger transition temperature of Tn/J = 0.590 
for the same systemii^ In Fig. ^ we present the magni- 
tude of the magnetization, m{T, B^), and the equilibrium 
angle, 9o{T, Bj.), as a function of the transverse field B^ 
for temperatures T below and above T/v. For all temper- 
atures, the magnetization m{T,Bx) increases when B^ 
is increased from zero. At the same time, 9o{T,Bx) also 
increases with increasing B^, indicating that the mag- 
netization is rotated towards the direction of the trans- 
verse field. Moreover, for T < T/v, the angle 6o{T,Bx) 
deviates from zero for infinitesimally small B^, implying 
that the rotation of the staggered moment does not re- 
quire a critical field strength, in contrast to the spin-flop 
transition associated with the application of a longitudi- 
nal magnetic field. '^■^■^ Note that for T > Tn, the hmit 
Br^ ^ leads to a non-zero angle 9o{T, B^) < 7r/2. This 
implies that a non-zero \mz(T, Bx)\ is induced by the 
transverse field, and that the component of the induced 
magnetization parallel to the transverse field increases 
faster than the component parallel to the easy axis. For 
T/J > 0.64 the behavior corresponds to the magnetiza- 
tion of an isotropic AFM.-^ 

The reorientation field Sfl(r) = 4m(T, Br){2J + D) is 
the smallest field at which the magnetization is parallel to 
the direction of the magnetic field, i.e., 9q{T, Bb) = tt/2. 
Both m(r, Bx) and 9q{T, B^) exhibit a discontinuous be- 
havior at the reorientation field. Specifically, m(r, B^) 
jumps at Bii(T) to a smaller value and increases with fur- 
ther increasing B^ > Bfj(T). At the same time, the angle 
9o{T,Bj;) jumps from 9 < 7r/2 to 7r/2 at B^. Presently 
we cannot judge whether this discontinuous behavior is 
'real' or an artefact of the approximations used above. 

The staggered magnetization component \mz{T, Bx)\ 
along the easy (z-) axis as computed from m(T, B^) and 
9o{T,Bx) is shown in Fig. |21as function of B^. Though 
this magnetization component is perpendicular to the ap- 
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FIG. 3: (a) Energy gap ei{Bx) and (b) spinwave stiffness 
gi{Bx), cf. Egs. 1121141 . of the magnon dispersion relation 
ei{Bx) as function of the transversal field Bx for T = 0. The 
insets show the behaviors of ei{Bx) and gi (B^) near the reori- 
entation field Br{0)/J ~ 4.0162. For the other denotations 
we refer to Fig.0 



plied field, we find that it also increases when is in- 
creased from zero. This behavior is quite unexpected 
since the rotation of m(T, Bx) into the direction of the 
transverse field should lead to a decrease in \mz{T, Bx)\- 
After passing through a maximum, \mz(T, Bx)\ discon- 
tinuously vanishes at the reorientation field Bfi{T). At 
the same time, the magnetization component mx{T, B^) 
along the magnetic field increases linearly with B^, and 
saturates for B^ > Bfi{T). As mentioned, for T > Tjv a 
non-zero B^ induces a finite \mz(T, Bx)\ which increases 
continuously from \mz\ = in the disordered (paramag- 
netic) state for Bx = 0, exhibits a maximum, and then 
vanishes at the reorientation field. Moreover, in the limit 
J = 0, i.e., for a purely Tsing'-like exchange D where 
the transversal magnetization components are missing, 
we obtain that \'mz{T, Bx)\ does not exhibit a maximum, 
but decreases monotonically. 

The increase of \mz\ for small B^ can be attributed 
to the suppression of quantum and thermal fluctuations, 
whose strength depends on the form of the magnon exci- 
tations spectrum. At T = 0, quantum fluctuations lead 
to a reduction of the staggered magnetization \mz\ from 
its saturation value given by S* = 1/2. In order to inves- 



tigate the strength of quantum fluctuations, we consider 
the magnon dispersion, £i_2(k) as a function of the trans- 
verse field Bx at T = 0. Near the center of the Brillouin 
zone at k = wc expand the dispersion to obtain 



ei,2(k) =e°2(S.) +31,2(5,) k2 



(13) 



where 2 is the energy gap ('mass') and gi^2 the spin- 
wave stiffness of the dispersion. For fields smaller than 
the reorientation field _Bij(0), we find 



D 



2J + D 



(4m)2(2J + D)2 -^2 



eliBx) = 
e°2{Bx) = {4m)^ D {2J + D) + B. 



-(14) 
(15) 



In Fig. 13 we present e^Bx) and gi{Bx) for the lower- 
energy magnon branch. Both quantities increase as 
Bx is increased from zero and exhibit a maximum at 
Bx/ J ~ 0.85 which coincides with the location of the 
maximum in Im^j. The increase of the excitation gap 
and of the spin stiffness reduce the strength of fluctua- 
tions, and are thus directly responsible for the increase 
in jm^j. A decreasing anisotropy D/J yields a smaller 
gap, which increases the strength of the fluctuations, and 
as a result, the maximum of \mz{0, Bx)\ becomes more 
pronounced. For Bx Bji{0) the dispersion softens^Si 
however, the gap retains a small but still finite value at 
the reorientation field Bji{0). At the same time gi{Bx) 
exhibits a pronounced spike. These features are respon- 
sible for the discontinuous behavior of the magnetization 
at the reorientation field discussed above. Note that for 
T > T/v and for a finite Bx the dispersions ei^2(T, Bx) 
do not become 'soft' but exhibit as expected a gap even 
when the magnetization is completely aligned with the 
transverse field, i.e., for 9o{T, Bx) = 7r/2. 

The increase of \mz{T, Bx)\ at small Bx is pronounced 
for systems where quantum fluctuations are most impor- 
tant, i.e., for a small spin quantum number S and for a 
low spatial dimensionality. Our calculations show that 
with increasing S the relative maximum of \mz{T, Bx)\ 
becomes smaller, and is not present at all for S = 00, i.e., 
for classical spins. 

In Fig. 01 we plot \mz{T, Bx)\ as a function of the 
temperature for various magnetic fields. We define the 
temperature at which \mz\ vanishes as the reorientation 
temperature Tr[Bx), with Tb{Bx = 0) = Tm- Note that 
already a weak magnetic field leads to a that is sig- 
nificantly larger than T/y. With further increasing Bx 
the reorientation temperature decreases and vanishes for 
Bx/J > 4.02. 

In the remainder of this Section we will briefiy men- 
tion related results that were obtained in other systems 
and by using different theoretical approaches. A max- 
imum of \mz{T, Bx)\ has recently been obtained for an 
anisotropic AFM Heisenberg chain»i^ Here \mz{T, Bx)\ 
increases from the paramagnetic state \mz\ = 0, since 
such a chain does not exhibit an ordered state, and cor- 
responds thus to temperatures T > T/y for the 2D AFM 
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FIG. 4: Staggered magnetization component \mz{T,Bx)\ 
along the easy axis as function of the temperature for dif- 
ferent . For the other denotations we refer to Fig. 

as investigated in this study. Solving the latter system 
with a single-spin mean-field approximation, no maxi- 
mum of \mz{T^ Bx)\ is obtained, since the transversal 
spin terms Sf SJ + Sf Sj of Eq. , which cause the ob- 
tained behavior, are neglected by this method. Never- 
theless, if these terms are taken into account to some 
extent, such as within a two-spin mean- field approxi- 
mation (Oguchi theorySi), the properties of \mz{T, B.j;)\ 
are qualitatively reproduced. Finally, a maximum of 
\mz{T,\Bz\, Bj;) \ is also obtained for an antiferromagnet- 
ically coupled Heisenberg spin pair, a system which can 
be solved exactly. In this case, for B^ — a finite mag- 
netic order is induced by a small staggered magnetic field 
\Bz\ along the z-axis. 

IV. CONCLUSION 

We have studied the magnetization of a 2D square- 
lattice anisotropic (XXZ-) AFM in a transverse mag- 
netic field. A many-body Green's function approach has 
been applied, which is known to yield a good descrip- 
tion of the magnetization for the case of a ferromag- 



netic monolayer jSi The fact that our theoretical approach 
yields values for the sublattice magnetization at B = 
and the Neel temperature that are similar to those ob- 
tained in quantum Monte Carlo calculations^^ supports 
the validity of the theoretical method also for the 2D 
AFM. 

We showed that the staggered magnetization 
\mz(T,Bx)\ along the easy-axis perpendicular to 
the field increases for small B^ and exhibits a maximum 
before vanishing at the reorientation field. For T > T/v, 
we demonstrate that the transverse field induces a non- 
zero magnetization |m^(T, i?a:)| which is perpendicular 
to the applied field. We argue that the increase of 
\mz{T,Bx)\ for small B^ arises from changes in the 
magnon excitation spectrum, which in turn leads to a 
suppression of thermal and quantum fluctuations. 

The described behavior of the staggered magnetiza- 
tion of a 2D AFM in a transverse field can possibly 
be observed by, e.g., x-ray magnetic linear dichroism 
(XMLD), since this method is sensitive to the magnitude 
of the magnetization components. An interesting ques- 
tion is whether a finite magnetization component along 
the easy-axis for temperatures slightly above T/v as in- 
duced by a transverse magnetic field can be measured, 
cf. Figs. EH 

Note that typical magnetic fields of a few Teslas yield 
Zeeman energies much smaller than the exchange, B^ ^ 
J. Hence, when such a transverse field is applied to the 
AFM, the canting angle will be small. In contrast, if the 
AFM is coupled to an ordered ferromagnet (FM), the 
intrinsic field due to the strong interlayer exchange cou- 
pling at the AFM/FM interface is considerably larger. 
The resulting angle 9o{T,Bx) could then be sufficiently 
large such that the results presented above are observ- 
able. A particular interest in such FM - AFM interfaces 
has revived lately in relation to the exchange bias effecti^ 

Useful discussions with K. D. Schotte are gratefully 
acknowledged. P. J. J. likes to thank the IPCMS in 
Strasbourg, France, for the hospitality, and the Deutsche 
Forschungsgemeinschaft, Sfb 290, for financial support. 
D.K.M acknowledges financial support from the Alexan- 
der von Humboldt Foundation. 



1 T. Tamaribuchi and M. Ishikawa, Phys. Rev. B 43, R1283 
(1991). 

^ V. Yu. Irkhin, A. A. Katanin, and M. I. Katsnelson, Phys. 

Rev. B 60, 1082 (1999). 
^ A. Cuccoli, T. Roscilde, V. Tognetti, R. Vaia, and P. Ver- 

rucchi, Phys. Rev. B 67, 104414 (2003). 
* M. E. Zhitomirsky and T. Nikuni, Phys. Rev. B 57, 5013 

(1998). 

^ M. E. Zhitomirsky and A.L. Chernyshev, Phys. Rev. Lett. 

82, 4536 (1999). 
^ D. Petitgrand, S. V. Maleyev, Ph. Bourges, and A.S. 

Ivanov, Phys. Rev. B 59, 1079 (1999); A. V. Syromyat- 



nikov and S. V. Maleyev, Phys. Rev. B 65, 012401 (2002). 
^ S. Jin, T. H. Tiefel, M. McCormack, R. A. Fastnacht, R. 

Ramesh, and L. H. Chen, Science 264, 413 (1994). 
* R. Melzi, S. Aldrovandi, F. Tedoldi, P. Carretta, P. Millet, 

and F. Mila, Phys. Rev. B 64, 024409 (2001). 
^ R. J. Birgeneau et al., Phys. Rev. B 38, 6614 (1988); E. 

Manousakis, Rev. Mod. Phys. 63, 1 (1991). 
° N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 

(1966). 

^ M. E. Lines, Phys. Rev. 133, A841 (1964); M. Bander 
and D. L. Mills, Phys. Rev. B 38, R12015 (1988); R. P. 
Erickson and D. L. Mills, Phys. Rev. B 43, R11527 (1991). 



6 



D. V. Dmitriev and V. Ya. Krivnov, Phys. Rev. B 70, 
144414 (2004). 

1^ M. Holtschneider, W. Selke, and R. Leidl, Phys. Rev. B 
72, 064443 (2005). 

B. B. Beard, R. J. Birgeneau, M. Greven, and U. J. Wiese, 
Phys. Rev. Lett. 80, 1742 (1998); J. K. Kim and M. Troyer, 
Phys. Rev. Lett. 80, 2705 (1998). 
1^ H. Q. Ding, J. Phys.: Condens. Matter 2, 7979 (1990); S. 
S. Aplesnin, Phys. Stat. Sol. B 207, 491 (1998). 
S. V. Tyabhkov, Methods in the quantum theory of mag- 
netism, Plenum Press, New York, 1967; W. Nolting, 
Quantentheorie des Magnetismus, vol.2, B. G. Teubner, 
Stuttgart, 1986. 

P. Frobrich and P. J. Kuntz, Europ. Phys. J. B 32, 445 
(2003). 

P. Frobrich, P. J. Jensen, and P. J. Kuntz, Europ. Phys. J. 
B 13, 477 (2000); P. Frobrich, P. J. Jensen, P. J. Kuntz, 
and A. Ecker, Europ. Phys. J. B 18, 579 (2000). 
P. Frobrich and P. J. Kuntz, J. Phys.: Condens. Matter 
17, 1167 (2005). 

S. Schwieger, J. Kienert, and W. Nohing, Phys. Rev. B 
71, 024428 (2005). 

M. G. Pini, P. Politi, R. L. Stamps, Phys. Rev. B 72, 



014454 (2005). 

S. V. Tyablikov, Ukr. Mat. Zh. 11, 287 (1959). 
23 P. Henelius, P. Frobrich, P. J. Kuntz, C. Timm, and P. J. 

Jensen, Phys. Rev. B 66, 094407 (2002). 
2* L Junger, D. Ihle, J. Richter, and A. Kliimper, Phys. Rev. 

B 70, 104419 (2004). 

D. L. MiUs, Phys. Rev. B 40, 11153 (1989). 
2® Usually the equilibrium angle 6o{T, Bx) should be deter- 
mined by minimizing the free energy F{T, 9) as obtained by 
the Green's function approach. Unfortunately, since within 
the Tyablikov decoupling F{T, 9) behaves unphysically at 
elevated temperatures,-^ this approach cannot be applied 
here. 

2^ T. Oguchi, Progr. Theor. Phys. 13, 148 (1955). 

2® G. van der Laan, Phys. Rev. B 57, 5250 (1998); J. Liining, 
F. Nolting, A. SchoU, H. Ohldag, J. W. Seo, J. Fompeyrine, 
J.-P. Locquet, and J. Stohr, Phys. Rev. B 67, 214433 
(2003). 

2^ For recent reviews see, e.g.: J. Nogues and L K. SchuUer, 
J. Magn. Magn. Mater. 192, 203 (1999); M. Kiwi, ibid. 
234/3, 584 (2001); R. L. Stamps, J. Phys. D: Appl. Phys. 
33, R247 (2000). 



